Integrated mapping and characterization of the gene underlying the okra leaf trait in Gossypium hirsutum L

Diverse leaf morphology has been observed among accessions of Gossypium hirsutum, including okra leaf, which has advantages and disadvantages in cotton production. The okra leaf locus has been mapped to chromosome 15 of the Dt subgenome, but the underlying gene has yet to be identified. In this study, we used a combination of targeted association analysis, F2 population-based fine mapping, and comparative sequencing of orthologues to identify a candidate gene underlying the okra leaf trait in G. hirsutum. The okra leaf gene identified, GhOKRA, encoded a homeodomain leucine-zipper class I protein, whose closely related genes in several other plant species have been shown to be involved in regulating leaf morphology. The transcript levels of GhOKRA in shoot apices were positively correlated with the phenotypic expression of the okra leaf trait. Of the multiple sequence variations observed in the coding region among GrOKRA of Gossypium raimondii and GhOKRA-Dt of normal and okra/superokra leaf G. hirsutum accessions, a non-synonymous substitution near the N terminus and the variable protein sequences at the C terminus may be related to the leaf shape difference. Our results suggest that both transcription and protein activity of GhOKRA may be involved in regulating leaf shape. Furthermore, we found that non-reciprocal homoeologous recombination, or gene conversion, may have played a role in the origin of the okra leaf allele. Our results provided tools for further investigating and understanding the fundamental biological processes that are responsible for the cotton leaf shape variation and will help in the design of cotton plants with an ideal leaf shape for enhanced cotton production.